Constraining the shape of the CMB: a Peak- by-Peak analysis. 
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The recent measurements of the power spectrum of Cosmic Microwave Background anisotropics 
are consistent with the simplest inflationary scenario and big bang nucleosynthesis constraints. 
However, these results rely on the assumption of a class of models based on primordial adiabatic 
perturbations, cold dark matter and a cosmological constant. In this paper we investigate the need 
for deviations from the A-CDM scenario by first characterizing the spectrum using a phenomeno- 
logical function in a 15 dimensional parameter space. Using a Monte Carlo Markov chain approach 
to Bayesian inference and a low curvature model template we then check for the presence of new 
physics and/or systematics in the CMB data. We find an almost perfect consistency between the 
phenomenological fits and the standard A-CDM models. The curvature of the secondary peaks is 
weakly constrained by the present data, but they are well located. The improved spectral resolution 
expected from future satellite experiments is warranted for a definitive test of the scenario. 



I. INTRODUCTION. 

The recent observations of the cosmic microwave 
background (CMB) anisotropics power spectrum 
(0,0, i, 0, i, i, mM) have pre- 

sented cosmologists with the possibility of studying 
the large scale properties of our universe with unprece- 
dented precision. As is well known (see e.g. [T^), 
the structure of the theoretical CMB spectrum, given 
mainly by the relative positions and amplitude of the 
so-called acoustic peaks, is sensitive to several cosmo- 
logical parameters. The existing CMB data sets are 
therefore being analyzed with increasing sophistica- 
tion (see 13j and 14] for important advancements) in 
an attempt to measure the undetermined cosmolo gica l 
quantities. The most recent analyses of this kind (EE|, 

m El, 113, 113, 113,124], 54J, in, [13, Mm) 

have revealed an outstanding agreement between the 
data and the inflationary predictions of a flat universe 
and of a primordial scale invariant spectrum of adia- 
batic density perturbations. Furthermore, the CMB 
constraint on the amount of matter density in baryons 
uoi) is now in very good agreement with the indepen- 
dent constraints from standard big bang nucleosyn- 
thesis (BBN) obtained from primordial deuterium (see 
e.g. [22]) and consistent within 2-(j with those 

derived from the combined analysis of ^He and ^ Li 
f p^). Finally, the detection of power around the ex- 
pected third peak, on arc-minutes scales, provides a 
new and independent evidence for the presence of non- 
baryonic dark matter ( 28]). 

The data therefore suggests that our present cosmo- 
logical model represents a beautiful and elegant theory 
able to explain most of the observations. 

However, the CMB result relies on the assumption 
of a particular class of models, based on adiabatic, 
passive and coherent (see [s^) primordial fluctuations, 
and cold dark matter. In the following we refer to this 
class of models as A-Cold Dark Matter (A-CDM). 

This weak point, shared by most of the current stud- 
ies, should not be overlooked: it has been recently 
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shown, for example, that the ver> 
sion of gravity waves (see e.g . |48lJ . [4yjj or isocur- 
vature modes (H^, |5l|, jsof ) into the analysis can 
completely erase most of the constraints derived from 
CMB alone. 

Furthermore, since even more exotic modifications 
like quintessence (|3^), topological defects (J3p,pl|), 
broken primordial scale invariances (H^, |34j . [35]), 
extra dimensions (^^) or unknown systematics (just 
to name a few) can be in principle considered, one 
should be extremely cautious in making any definitive 
conclusion from the present CMB observations. 

It is therefore timely to investigate if the present 
CMB data are in complete agreement with the A- 
CDM scenario or if we are losing relevant scientific 
informations by restricting the current analysis to a 
subset of models (see e.g. |3^). 

In the present paper we check to what extent modifi- 
cations to the standard A-CDM scenario are needed by 
current CMB observations with two complementary 
approaches: First, we provide a model-independent 
analysis by fitting the data with a phenomenologi- 
cal function and characterizing the observed multi- 
ple peaks. Phenomenological fits have been exten- 
sively used in t he p ast and recent CMB analyses 
(in, IHi, 0, lU, EH, lEl). Our analysis dif- 
fers in two ways: we include the latest CMB data 
from the Boomerang (ll^ ) , VSAE ([TJ), ACBAR 
(1^), and Archeops (I§|) experiments and we make 
use of a Monte Carlo Markov Chain (MCMC) algo- 
rithm, which allows us to investigate a large number 
of parameter simultaneously (15 in our case). 

We then compare the position, relative amplitude 
and width of the peaks with the same features ex- 
pected in a 4-parameters model template of A-CDM 
spectra. By doing a peak-by-peak comparison be- 
tween the theory and the phenomenological fit which 
is based on a much wider set of parameters, we then 
verify in a systematic way the agreement with the 
standard theoretical expectations. 

As a by-product of the analysis, we present a set 
of cosmological diagrams that directly translate, un- 
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der the assumption of A-CDM, the constraints on the 
features in the spectrum into bounds on several cos- 
mological parameters. These diagrams offer the op- 
portunity of quick, by-eye, data to model comparison. 

Our paper is organized as follows: In section II 
we discuss the phenomenological representation of the 
power spectrum, the analysis method we used and the 
MCMC algorithm. In section III we present our re- 
sults. Finally, in section IV, we discuss our conclu- 
sions. 



II. PHENOMENOLOGICAL 
REPRESENTATION. 



We model the multiple peaks in the CMB angular 
spectrum by the following function: 



N 



i{i + l)C£/27r = AT^ exp(- 



■QV^^^) (1) 



where, in our case, N = 5. In order to avoid de- 
generacy of overlapping gaussians, we parametrise the 
centers of the secondary gaussians as functions of the 
positions of the previous gaussians: 

h = h{l^a) (2) 

Is = + a + 

h = /i(l + a + /3 + 7) 

k = + a + /3 + 7 + (5) 



We use this formula to make a phenomenological fit 
to the current CMB data, constraining the values of 
the 15 parameters ATj, ii and cr^. 

The use of gaussian-shaped function to describe the 
CMB spectrum is now becoming a standard method 
in the literature (see e.g. 39], 41], 0, 0). A ma- 
jor difference with respect to previous works is that 
we are using only one fitting function, varying all its 
parameters simultaneously, while in general peaks are 
characterized with one single function in different se- 
lected regions of the spectrum, in correspondence with 
the expected peaks. 

The advantage of a single fitting function is a better 
control of the correlations between the phenomenolog- 
ical parameters as we show in the next section where 
we report the values of the correlation matrix. 

Recently, Douspis and Ferreira ( 52]) used a Gaus- 
sian plus an oscillating function as a phenomenologi- 
cal model. The method used here is more general, in 
the sense that we allow independent amplitudes and 
widths of the secondary peaks as well as impose no 
periodicity. 

We use the CMB data as listed in tabled spanning 
the range 50 < £ < 2500. 



For all the experiments we use the publicly available 
window functions and correlations in order to com- 
pute the expected theoretical signal Cb inside the bin. 
The likelihood for a given phenomenological model is 
defined by -2\nC = {Cf - Cf)MBB'{Cf, - Cf>) 
where Mbb' is the Gaussian curvature of the likeli- 
hood matrix at the peak. When available, we use the 
lognormal approximation to the band-powers. 

We marginalize over the reported Gaussian dis- 
tributed calibration error for each experiment and 
we include the beam uncertainties by the analytical 
marginalization method presented in ("S^). 

We perform our analysis on 3 different data sets: 
the full data set, the low- frequency (LF) data (experi- 
ments that covered frequencies in the electromagnetic 
spectrum of 90 GHz and below) and a high-frequency 
(HF) data (frequencies higher than 90 GHz). The 
reasons of this choice are twofold: First, any discrep- 
ancy between the 2 analyses would hint towards the 
presence of undetermined foreground. Secondly, this 
facilitates the comparison with future and soon to be 
released observed power spectra at "low" frequencies, 
like those expected from the MAP satellite. We test 
the stability of our result by including a set of older 
CMB experiments as reported in tableland we also 
check that our model contains no bias towards the 
presence of peaks by fitting a set of mock data from a 
spectrum (see fig. EJ. 
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FIG. 1: The CMB data used in this analysis: beam and 
calibration errors are not included. 



The phenomenological fit is operated through an 
MCMC algorithm. The MCMC approach is to gen- 
erate a random walk through parameter space that 
converges towards the most likely value of the param- 
eters and samples the parameter space following the 
posterior probability distribution. In the general case, 
the number of parameters and their priors have to be 
defined. There is no limit in resolution (except nu- 
merical precision of the computer). The priors define 
the volume in parameter space in which the random 
walk takes place. 

At each iteration, a point x^+i is randomly selected 
in the m-dimensional parameter space. Its likelihood 
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TABLE I: List of experimental data used in this study. We examine low and high frequency experiments together and 
separately. We then check our results by including some older experiments. 



Experiment 




/ range 


reference 


High Frequency Experiments (HF): 


Acbar 




150 - 


- 3000 


Kuo et a/., 2002 


Archeops 




15 - 


- 350 


Benoit et a/., 2002 


Boomerang 98 




50 - 


1000 


Ruhl et al, 2002 


Maxima 




73 - 


1161 


Lee et al., 2001, Ap. J. 561 (2001) L1-L6 


Low Frequency Experiments 


(LF): 








DASI 




117 


- 836 


Halverson et a/., Ap. J, 568, 38, 2002 


CBI 




400 - 


- 1450 


Pearson et a/., 2002 


VSAE 




160 - 


- 1400 


Grainge et a/., 2002 


Older Experiments: 


TOCO 97 




63 - 


- 194 


Miller et a/.Ap. J. Supp., 140, 115, 2002 


TOCO 98 




128 


-409 


Miller et a/.Ap. J. Supp., 140, 115, 2002 


MSAM 




90 - 


-400 


Wilson, et al., 2002 


QMASK 




105 


- 359 


Xu et al, Phys. Rev. D65, 2002 


Python V 




67- 


- 267 


Coble et al, 2001 


Boomerang NA 




50 - 


-400 


P. Mauskopf et a/., Ap. J. Letters, 536, L59, 2000 



the best fit value is obtained by averaging over the 
samples. The MCMC procedure is described in more 
detail in IH and l^. 

We use wide uniform priors for each parameter, in 
which our resulting likelihoods are fully encompassed, 
except for the widths of the secondary gaussians which 
are not well constrained by the data. We run the 
MCMC routine in order to get ~ 7000 samples for 
each subset of data. Also, in order to check that out 
parametrisation is not biased towards the presence of 
peaks, we test it to fit a flat spectrum. We generate 
a set of mock data by convolving a power spectrum 
consisting of a peak and a flat line with the exper- 
imental window functions. The results show that a 
flat power spectrum is easily recovered within our pri- 
ors, as shown in figure [21 

We then consider a flat, adiabatic, A-CDM model 
template of CMB angular power spectra, computed 
with CMBFAST sampling the various param- 

eters as follows: Qcdmh'^ = ^cdm = 0.01, ...0.40, in 
steps of 0.01; Vtbh^ = ujj, = 0.001, 0.040, in steps 
of 0.001 and = 0.0, ...,0.95, in steps of 0.05. The 
value of the Hubble constant is not an independent 
parameter, since: 

We vary the spectral index of the primordial density 
perturbations within the range Us = 0.60, ...,1.40 (in 
steps of 0.02). 

For each model in the template we then consider the 
corresponding values AT^, such that the formula 
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FIG. 2: The 10 best MCMC samples fitting mock data. 
The data were generated by convolving the experimental 
window functions with a power spectrum consisting of a 
first peak followed by a flat line. 



is evaluated by comparing to the data. A sample is 
said to be accepted into the Markov chain or not, de- 
pending on the following acceptance criterion: 



u < min 



injl, 



p(x„+i)£(Xn+i) 

p(x„)£(x„) 



(3) 



where is a random number sampled from a uniform 
distribution on the interval [0,1], p(x) is the prior 
and jC{x.) is the likelihood, containing the information 
from the data. At the end of the MCMC routine, the 
samples are counted and their number density pro- 
jected onto one or more dimensions is proportional to 
the marginalised posterior distribution of the param- 
eters. If the posterior follows a gaussian distribution. 
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in Eq. (1) represents the best fit to its shape. Indeed, 
the ATi do not exactly correspond to the amplitudes 
of the peaks, as our spectrum is a sum of gaussians 
and the power from each gaussian contribute over the 
whole spectrum. We find that, restricting the range in 
^ to 50, 1500, equation (1) approximates the shape 
of the spectra in our template well (better than ~ 10% 
in Ce). 

We also check that the use of different phenomeno- 
logical functions such as lorentzians or log-normals has 
no relevant effect on our results. 

It is important to note that we restrict the 
parametrisation of our template of theoretical mod- 
els to a set of only 4 parameters. However, because of 
the 'cosmic degeneracy' in the CMB observables, this 
is enough to describe the possible shapes of the CMB 
spectra in the A-CDM scenario. Increasing the optical 
depth Tc or adding a background of gravity waves, for 
example, is nearly equivalent to changing some of the 
parameters already considered like ns and uj^. On the 
other hand we characterize the peaks in the spectrum 
with a phenomenological fit based on 15 parameters, 
which allows independent positions, amplitudes and 
widths of the observed features. 

The comparison between the model-independent 
values ATi, and cr^ obtained by fitting the data with 
Eq. (1) and the corresponding values expected in the 
template of theoretical models represents therefore a 
strong check of the theory and can give hints for the 
presence of systematics and/or new physics. 



III. RESULTS 

In Table HH we report the 68% limits on AT^, £i 
and ai of Eq. 1 obtained by analyzing the present 
CMB data with an MCMC procedure for each sub- 
set of CMB data. We also report the constraints on 
combinations of those parameters that are more di- 
rectly connected to the cosmological parameters (see 
the discussion below). 

In Table IIIII we report the correlation matrix be- 
tween the parameters of our phenomenological fit. As 
we can see, important correlations exist between the 
parameters: for example, the amplitude of the peaks 
is highly anti-correlated with the widht of the adia- 
cent gaussians. This further illustrate the utility of 
analyzing the data with a single fitting function in or- 
der to properly evaluate the statistical significance of 
the oscillations. 

Figure |31 shows the marginalised likelihood func- 
tions of the important CMB observables, and the 
agreement between the three data subsets. Figure 0] 
shows the best-fit cosmological and phenomenological 
models. Although our sums of gaussians comprises 
any cosmological power spectrum to within 10%, their 
preferred shape differs from the cosmological model 
owing to the large parameter space allowing to fit the 
data very well. 



Figure [51 shows 68% and 95% confidence levels in 
power in the i-AT'^ plane. This shows the most likely 
path of the power spectrum. These figures are gen- 
erated as follows: First, we remind the reader that 
the number density of MCMC samples is proportional 
to the marginalised posterior. We divide the i-AT'^ 
plane into pixels of size 10 x 10. For each pixel, we 
count the phenomenological spectra that go through 
it. Hence, this is a density plot of our MCMC sam- 
ples. The absence of a 68% region at certain i ranges 
indicates that more data are needed to constrain the 
power. 

From figure [5J it appears that LF experiments pro- 
vide constraints at high multipoles better than HF 
experiments which constrain the power at intermedi- 
ate scales strongly. HF experiments also provide tight 
constraints at low ^, mainly due to Archeops. 

We also show the same result obtained by using the 
old data only. We find that the analysis provides no 
constraints above £ ~ 400, as expected from the data. 
It is also consistent with the Archeops observations. 
Hence, when combining the old data with the LF data, 
HF data or both, our results remain essentially the 
same. 

The values are in reasonable agreement with the 
results obtained by similar analyses (see e.g. 0, [Zlj . 
[52I, [22) and point towards the presence of multiple 
peaks in the spectrum. 

The low frequency and high frequency experiments 
yield consistent results, showing that possible system- 
atics due to galactic foregrounds are under control. 
However it is worth noticing that the LF data are more 
consistent with higher amplitude of secondary peaks. 
These experiments also use different techniques. The 
HF experiments are bolometers whereas the LF ex- 
periments are interferometers. The different nature of 
experimental uncertainties as well as their evaluation 
might suffer from different weaknesses contributing to 
enhance this contrast. 

The CBI and ACBAR data at high i are in agree- 
ment with the expected damping tail (see e.g. 61]). 
However, the poor spectral resolution (A^ ~ 200) does 
not allow us to constrain subsequent peaks. 

Still, it is interesting to compare the values obtained 
with those expected in the A-CDM scenario with dif- 
ferent priors, as we do in the last three columns of 
Table m (see caption). 

Generally speaking, considering that the reported 
errors are at 1-a and that the theoretical models are 
COBE normalized, which allows a further 10% shift in 
amplitude, the model-independent values are in very 
good agreement with those predicted by the A-CDM 
models. It is interesting to notice that for the full data 
set, the subsequent peaks appear to be slightly lower 
in amplitude than those expected in the concordance 
model with = 1. This favors a spectral index ns 
slightly lower than one. However, there is a strong 
degeneracy between ns and the optical depth Tc (see 
e.g. ^^) and models with = 1 can be put in better 
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agreement with the observations when increasing Tc- 
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FIG. 3: Marginalised likelihood functions for CMB observ- 
ables. The red (solid), green (dashed) and blue (dotted) 
curves correspond to the following three data sets, respec- 
tively: Fuh data set, LP, HF. 
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FIG. 4: Comparison of the concordance cosmological 
power spectrum and the best-fit phenomenological power 
spectrum for each data subset.). 




FIG. 5: 68% and 95% confidence levels in the AT^, ^ plane. 
The first three panels correspond to the full data set, HF 
data and LF data, respectively. The last panel shows the 
same result obtained by analysing the old data only. It is 
found to be consistent with the more recent data, and pro- 
vides no constraints above I ^ 400 where the data provide 
no measurements. This shows where the power is well con- 
strained by the data using our phenomenological function. 
Where no limits of power are found within 68%, it means 
that the data are not sufficient to constrain the power at 
that angular scale. 



We further investigate the consistency with A-CDM 
by considering phenomenological diagrams relating 
the relative amplitudes and positions of the peaks with 
variations in a specific physical parameter. 

In the A-CDM adiabatic scenario, two key param- 
eters control the relative power between the first and 
second peaks: the physical baryon density cOb and the 
primordial spectral index (see e.g. jH^I)- Increasing 
uji) enhances the odd-numbered peaks relative to the 
even-numbered ones, while increasing ns enhances the 
small-scale peaks relative to the ones at larger scales. 



In figure we plot the values allowed in our model 
template (restricted by a set of rather conservative 
cosmological constraints, see the caption) for the rel- 
ative amplitude AT1/AT2 as functions of the param- 
eters ujj) and ns. As expected, increasing (decreas- 
ing) LOi) (ns) increases AT1/AT2. The region is very 
broad, mainly owing to the degeneracy between these 
two parameters. Nevertheless, super-imposing the 1-cr 
phenomenological constraint on AT1/AT2 in the dia- 
gram provides interesting constraints. Models with 
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TABLE II: First 3 columns: 1-a constraints on the parameters of the phenomenological model for three subsets of 
data: The full data set, LF data and HF data. The allowed range for the same parameters in a database of COBE 
normalized theoretical models (no experimental data are considered) is also reported in the following 2 columns for the 
case of weak priors (0.05 < Qcdm < 0.5, 0.15 < Qbh'^ < 0.25, 0.55 < h < 0.88, 0.80 < ns < 1.10) and strong priors 
(0.10 < Qcdm < 0.35, 0.18 < Qbh^ < 0.22, 0.65 < h < 0.80, 0.95 < < 1.05). In the last column we also show the 
values for a COBE normalized concordance model with Qcdm = 0.31, Qb = 0.04, h = 0.7 and Us = 1. 
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Full data set 


High Frequency 


Low Frequency 


A-CDM 


A-CDM 


Concordance 


Observable 




Experiments 


Experiments 












Phen. Fit 


Phen. Fit 


Phen. Fit 


Weak Priors 


Stron 


g Priors 




ATi 








(48- 110)/iK 


(65- 


95)/iK 


74.9/xK 


AT2 








(32 - 76)i^K 


(46 - 




b3.QpK 


ATs 






54Atlt.2f^K 


(36 - 67)i^K 


(49- 


m)yiK 




AT4 






27.9l?;^M^ 


{22 - 40)fiK 


(26 - 


37)fiK 


35.6fjK 


AT5 




18 Atlii-iK 




(19 - 36)^K 


(22 - 


?>2)iiK 


27.6//K 


ii 


208.81^;? 


204.61^^9^ 


206.8l^°:g 


192 - 279 


203 


-242 


215 


i2 






rqq+25 
000_2Q 


464 - 696 


488 - 594 


514 


is 








692 - 1086 


732 


-919 


781 


£4 




1158+^f 


11891^7 


1140 - 1386 


1210 


- 1301 


1190 


£5 


1474l^f 




I515I34Q 


1380 - 1590 


1460 


- 1550 


1491 


CTl 


93.3lt.2 




QQ + 12.7 


86 - 136 


86- 


- 108 


93 


(72 


111 2+^'^ '^ 

' unbounded 


78.2l}^.3 


ai Q+36.7 

^unbounded 


71 - 121 


71 - 


- 107 


86 


crs 


09 ^+20. 7 


J.OU.O„„|,ounded 


Ut^.O^nbounded 


86 - 150 


86- 


- 150 


102 


(74 


not constrained 


not constrained 


not constrained 


70 - 120 


70 - 


- 110 


87 


(75 


not constrained 


not constrained 


not constrained 


70 - 120 


75- 


- 110 


88 


h/i2 


n QQ9+0.033 
U.OOZ_Q 


n 415+0014 


U.O/ U_Q 


0.379 - 0.429 


0.402 


- 0.422 


0.418 


h/is 


U.ZOO_Q q;|^q 






0.242 - 0.290 


0.259 


- 0.281 


0.275 


AT1/AT2 


-L.OO_Q Q7 




l-44l°:l«5 


1.23 - 1.68 


1.32 


- 1.51 


1.40 


ATi/ATs 


-'-•'-'^-0.11 


-L-00_Q j^O 


1 2q+°-22 


1.12-2.13 


1.26 


- 1.73 


1.37 


AT2/AT3 


-'-•^-'--0.09 


QQ+O-20 




0.80 - 1.32 


0.92 


- 1.17 


0.98 


i2-ii 




290t?9 


329ttl 


267 - 427 


284 


-353 


299 


is-ii 


611139 


5551^1 




484 - 817 


528 


-679 


566 



TABLE III: Correlation matrix for the 15 phenomenological parameters when fitted to the full data set. 





ATi 


AT2 


AT3 


AT4 


AT5 


ii 


i2 


is 




4 


(71 


(72 


(7s 


(74 


(75 


ATi 


1.00 


0.20 


0.32 


0.12 


0.06 


-0.02 


-0.17 


-0.07 


0.03 


-0.02 


-0.23 


-0.19 


0.12 


-0.18 


0.04 


Ar2 


0.20 


1.00 


0.03 


-0.06 


-0.07 


-0.30 


0.54 


0.54 


-0.15 


0.00 


-0.32 


0.38 


-0.70 


0.18 


0.11 


ATs 


0.32 


0.03 


1.00 


0.37 


0.07 


0.38 


-0.54 


-0.27 


0.40 


0.08 


0.30 


-0.63 


0.40 


-0.45 


-0.08 


AT4 


0.12 


-0.06 


0.37 


1.00 


0.32 


0.06 


-0.19 


0.08 


0.20 


-0.10 


0.03 


-0.13 


0.27 


-0.71 


-0.10 


AT5 


0.06 


-0.07 


0.07 


0.32 


1.00 


-0.03 


0.01 


0.16 


-0.36 


-0.61 


-0.07 


0.05 


0.06 


-0.51 


-0.09 


ii 


-0.02 


-0.30 


0.38 


0.06 


-0.03 


1.00 


-0.42 


-0.55 


0.10 


0.08 


0.60 


-0.75 


0.49 


-0.09 


-0.05 


i2 


-0.17 


0.54 


-0.54 


-0.19 


0.01 


-0.42 


1.00 


0.72 


-0.37 


-0.09 


-0.32 


0.79 


-0.84 


0.31 


0.08 




-0.07 


0.54 


-0.27 


0.08 


0.16 


-0.55 


0.72 


1.00 


-0.04 


-0.13 


-0.52 


0.77 


-0.58 


-0.07 


0.01 


i4 


0.03 


-0.15 


0.40 


0.20 


-0.36 


0.10 


-0.37 


-0.04 


1.00 


0.64 


0.10 


-0.23 


0.48 


-0.20 


-0.25 




-0.02 


0.00 


0.08 


-0.10 


-0.61 


0.08 


-0.09 


-0.13 


0.64 


1.00 


0.11 


-0.09 


0.10 


0.34 
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FIG. 6: Observed relative amplitude between the first and 
second peak and comparison with adiabatic A-CDM cos- 
mological models. The red region ("A-CDM no CMB") 
defines all the values of AT1/AT2 present in the model 
template described in the text with the additional priors 
0.1 < Qcdm < 0.5 and 0.55 < h < 0.88. In the top panel, 
the additional prior 0.015 < Qfe/i^ < 0.025 is included, 
while in the bottom panel we use 0.8 < ns < l-l- The 
68% c.l. constraint obtained by the phenomenological fit 
and he 95% c.l. of all the A — CDM models compati- 
ble with CMB are also plotted for comparison. The con- 
straints on Qbh^ from standard BEN from D (Buries et 
al.2001) and ^He +^ Li (Cyburt et al. 2001) observations 
are also included in the second plot. 



a value of the spectral index > 1.15, which are 
not easily accommodated in most inflationary models 
(see e.g. and references therein) or in disagree- 
ment with the BEN constraint Ub = 0.020 ± 0.002 
are in fact not favored by our model-independent fit. 
In that figure, we also plot the constraints obtained 
by fitting the CME data with the models in the tem- 
plate. This reduces in a severe way the number of 
allowed A-CDM models. Nevertheless, the result on 
relative amplitude of the peaks is completely consis- 
tent with the one derived by the phenomenological 
fit. This method provides better constraints on the 
amount of cold dark matter cocdm if we consider the 
relative amplitudes AT1/AT3 and AT2/AT3 as we do 
in the top and bottom panels of figure [7| A decrease 
in ujcdm has the effect of decreasing the amplitude of 




0.6 - 
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FIG. 7: Observed relative amplitudes between the first and 
third peaks, and between the second and third peaks, and 
comparison with adiabatic A-CDM cosmological models. 
The red region ("A-CDM no CMB") defines ah the values 
of AT1/AT3 present in the template of models described 
in the text with the additional priors 0.55 < h < 0.88, 
0.015 < Qb/i^ < 0.025 and 0.8 < ns < 1.1. The 68% c.l. 
constraint obtained by the phenomenological fit and the 
95% c.l. of ah the A-CDM models compatible with CMB 
are also plotted for comparison. 



the third peak (see e.g. 57]). As we can see, the two 
observational values of AT1/AT3 and AT2/AT3 pro- 
vide similar constraints on cOcdm in a non trivial way 
with iVcdm ^ 0.12 and with cVcdm = or Ucdm > 0.3 
in disagreement with the data. 

In a flat A-CDM model a variation in Qa shifts the 
spectrum as ^ — > IZi with the shift parameter IZ given 
by (see 



rZdec 

n=^/\^Q [{i-nA){i + zf + nA]-'^^dz. 

Jo 

(5) 

Varying the Hubble constant, parametrized as Hq = 
lOOh Km/sec/Mpc, changes the scale of equality and 
produces a similar shift. These two parameters are 
related to the age of the universe by: 



8 



to = -9.SGy / [h\{l - (^a)(1 + zf + QaT^^^cIz 

Jini 

(6) 

We can therefore expect that a deterniination of 
the peak positions provides constraints on these three 
quantities. 

In figure IHl and figure El we plot similar diagrams as 
above for £2 — ii and is — ii as functions of (^a, h and 
age, to. Even if the region of the allowed models is 
quite broad because of the intrinsic degeneracies, the 
observed peak positions strongly favor a model with 
cosmological constant > 0, a Hubble parameter 
h < 0.8, compatible with the Hubble Space Telescope 
(HST) result of h = 0.72±0.08 dHU), and an age to > 
14 Gy, compatible with the age of the oldest globular 
clusters (see e.g. 6OJ). Again, the values obtained by 
the phenomenological fit are in agreement with those 
derived by the standard CMB+A-CDM analysis. 

Another parameter that affects the position of the 
peaks is the spectral index n^. However, the effect 
is different, the shift being scale dependent. There- 
fore, it is better to consider the quantities ^1/^2 and 
ii/is which are unaffected by the overall shift IZ. The 
corresponding diagrams are plotted in figure 11171 As 
mentioned earlier, the observed values point towards 
a low value of the spectral index < 1. 



IV. CONCLUSIONS 

In this paper we investigated the consistency of the 
most recent CMB data with a class of A-CDM adia- 
batic infiationary models. First we characterized the 
positions, amplitudes and widths of the peaks by fit- 
ting the data with simple phenomenological functions 
composed by several gaussians. The detection of the 
peak amplitudes and positions is quite robust and sta- 
ble between different data sets. We found that all 
the features are consistent with those expected by the 
standard theory. We also examined where the data 
contains the most information in the power-angular 
scale plane. We found that the low frequency ex- 
periments provide good constraints at small angu- 
lar scales, consistent with the expected damping tail, 
whereas high frequency experiments provide strong 
limits on the power at large and intermediate scales. 
We observe that HF experiments and LF experiments 
yield very consistent results, although LF data seems 
to provide evidence for higher secondary oscillations. 
Overall, the power spectrum is now well determined 
until £ ~ 1500. The inclusion of older data does not af- 



fect our conclusions as they do not measure the power 
beyond ^ - 400. 

Furthermore, we related the features in the spec- 
trum with several cosmological parameters by intro- 
ducing cosmological diagrams that can be used for 
quick, by-eye, parameter estimations. 

The relative amplitude of the first and second peak, 
in particular, of about ~ 1.56 is consistent with the 
baryon density expected from BBN and suggests a 
value of lower than one in the case of negligible 
reionization. The amplitudes of the third peak rela- 
tive to the first and to the second, AT1/AT3 ~ 1.6 
and AT2/AT3 ~ 1 strongly suggest the presence of 
cold dark matter but also limits time its contribution 
to values Ucdm < 0.2. The relative positions of the 
peaks, £2 — £-1 ^ 330 and £3 — ii ~ 610 is point- 
ing towards the presence of a cosmological constant, a 
Hubble parameter on the low side of the value allowed 
by the recent HST measurements {h ~ 0.65) and to 
an age of the universe to ~ 14.5 Gyrs consistent with 
the measurements of the oldest globular clusters. 

It is reassuring that all those conclusions, obtained 
by just drawing few lines in the diagrams presented 
in Figs. 5 — 9, are in agreement with the results ob- 
tained by a more careful standard analysis. Within 
the models considered in our database we found (at 
68% C.I.): Us = 0.96±0.03, ujb = 0.022±0.003, ujcdm = 
0.12±0.03, = 0.63±0.16, and to = 14.2±0.7 Gyrs. 

The results obtained here show no need for 
modifications to the standard model, like gravity 
waves, quintessence, isocurvature modes, or extra- 
backgrounds of relativist ic particles. Furthermore, 
possible systematic effects due to unknown fore- 
grounds or calibration and beam uncertainties are not 
immediately suggested, since the different data sets 
are consistent with the theory. 

Even if the width of the gaussians is poorly con- 
strained, we found supporting evidence for multiple 
oscillations in the data between 430 < £ < 910. Be- 
yond that, the newest experimental results show a 
damping of the power. 

It is the duty of future satellite CMB experiments 
to point out discrepancies that might place the possi- 
bility of new physics in a more favorable light. 
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Age (Gyr) 

FIG. 8: Observed relative positions between the first and 
second peak and comparison with adiabatic A-CDM cos- 
mological models. The red region ("A-CDM no CMB") 
defines all the values of ^2 — ii present in the model 
template described in the text with the additional pri- 
ors 0.015 < Qbh'^ < 0.025 and 0.8 < ns < 1.1. In 
the top panel (constraints on Qa) the additional prior 
0.55 < h < 0.88 is used. In the centre panel, we use 
0.1 < ^cdm < 0.5. The band on the ^-axis on the bot- 
tom panel is the constraint on the age of the oldest halo 
globular cluster in the sample of Salaris and Weiss (1998). 
The 68% c.l. constraint obtained by the phenomenological 
fit and the 95% c.l. of all the A-CDM models compatible 
with CMB are also plotted for comparison. 
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FIG. 9: Observed relative positions between the first and 
third peak and comparison with adiabatic A-CDM cosmo- 
logical models. The definitions of the regions and priors 
used are the same as in Fig. El 
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FIG. 10: Observed relative positions between the first, 
second and third peak and comparison with adiabatic A- 
CDM cosmological models. The red region ("A-CDM no 
CMB") defines all the values of ^1/^2 and ^1/^3 present 
in the modeltemplate described in the text with the ad- 
ditional priors 0.015 < Qbh'^ < 0.025, 0.55 < h < 0.88 
and 0.1 < ^cdm < 0.5. The 68% c.l. constraint obtained 
by the phenomenological fit and the 95% c.l. of all the 
A — CDM models compatible with CMB are also plotted 
for comparison. 
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